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The Feynman propagator used in the conventional in-out formahsm in quantum field the- 
ory is not a causal propagator as wave packets are propagated virtually instantaneously 
outside the causal region of the initial state. We formulate a causal in-out formalism in 
quantum field theory by making use of the Wheeler propagator, the time ordered commu- 
tator propagator, which is manifestly causal. Only free scalar field theories and their first 
quantization are considered. We identify the real Klein Gordon field itself as the wave func- 
tion of a neutral spinless relativistic particle. Furthermore, we derive a probability density 
for our relativistic wave packet using the inner product between states that live on a suit- 
ably defined Hilbert space of real quantum fields. We show that the time evolution of our 
probability density is governed by the Wheeler propagator, such that it behaves causally 
too. 

PACS numbers: 03.70.+k, ll.lO.-z, 03.65.Pm, 11.40.-q, 11.80.-m, 31.30.jd 

I. INTRODUCTION 

Quantum field theory is a very successful theory. The Feynman rules are very simple to follow 
but stem from complicated underlying considerations. For example, in order to compute a certain 
transition amplitude, a particular Feynman diagram needs to be renormalized and computed before 
it can be related to observables (a cross section). The standard in-out formalism used to calculate 
Feynman diagrams has however two shortcomings: 

• Transition amplitudes in the in-out formalism are not causal; 

• Transition amplitudes in the in-out formalism are not real. 
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The standard in-out formalism is not causal in the following sense: the Feynman propagator that 
one uses to calculate diagrams in perturbation theory is not a causal propagator as it has non- 
vanishing support outside the past- and future light cone of a certain spacetime point x. In other 
words, the Feynman propagator /S.f{.x;x') does not vanish when x and x' are spacelike separated. 
This sense of causality is sometimes referred to as Einstein causality. 

The standard in-out formalism is not real in the following sense: even for real scalar fields 
the quantum description is complex since creation operators only create the positive frequency 
part of a field and all negative frequency dependence is automatically projected out. Self-energies 
are therefore generically complex valued although we started the computation with real fields. In 
other words, the in-out effective action does not generate the correct quantum corrected equations 
of motion for the real scalar field 

HQ. 

Similar to unitarity Q], causality is an ingredient of fundamental importance of relativistic 
quantum field theories. Let us be very explicit about the following point: no one would claim that 
quantum field theory is not a causal theory. Expectation values can be expressed in terms of series 
of time ordered nested commutators (see e.g. P, Q]), which makes them manifestly causal. We 
thus conclude that when calculating observables using Feynman propagators crucial cancellations 
must occur between various diagrams such that the final result is causal. 

The approach which computes true expectation values without computing scattering amplitudes 



as an intermediate step is referred to as the in-in or Schwinger-Keldysh formalism [^, I6l4l0l|. for 



applications see 
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IJ]. The in-in formalism does not suffer from the shortcomings of the in-out 
formalism as mentioned above^. Also, self-mass corrections to real propagators are purely real as 
they should. Hence, the interpretation of the self-mass corrections thus obtained is much more 
straightforward than in the in-out formalism. 

The main idea that we outline in this paper is how to set up a real and causal first quantized 
description in perturbative quantum field theory in an in-out setting. We only develop the ingre- 
dients necessary to build a perturbative quantum field theory that meets these requirements, and 
we postpone the development of a perturbation theory itself to a future paper. Let us begin by 
recalling the first quantized picture in conventional quantum field theory before we explain where 
we differ from the standard approach. Quantum states in the conventional approach to quantum 
field theory are complex valued, just as in non-relativistic quantum mechanics. The reason is very 



^ Although the in-in formalism is perfectly causal, it is not manifest in all incarnations of the in-in formalism. The 
diagrammatic approaches of [l^ and [Tg] are for instance not based on purely causal propagators. 
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simple to see in second quantization: a creation operator creates a positive frequency wave mode 
out of the vacuum ajj|0) = \ijJk,k), where ijJk = ij? + rm?)^/"^ as usual. Now one can employ the 
superposition principle such that the positive frequency contribution 4>-\-{x) to the real field 4>{x) 
is a quantum mechanical wave function: 

(1) 

Clearly, this wave function is complex valued. Furthermore, by making use of the inner product 
between fields that live on the Hilbert space of positive frequency solutions to the Klein Gordon 
equation one can straightforwardly find a probability density that, integrated over its domain of 
non-zero support, can be properly normalized to unity. See for example section 3.1 of {l?! or section 
14.2 of [l^ where a similar interpretation is made. This probability density is the 0th component 
of the familiar Klein Gordon current given by: 

iVo^^(x) = <^;(x)i5t</.+ (x), (2) 

The natural question that presents itself when writing down an equation like ([T|) and ([2]) is: why 
is ^(x) not the wave function in quantum field theory? In other words: why is ^{x) not the 
relativistic generalization of the wave function of ordinary quantum mechanics? The answer is 
that the candidate wave function ^{x) defined above does not behave causally (see e.g. [19l . 



[20I). Advanced or retarded propagation of the candidate wave function ^'(x) is governed by the 
Feynman propagator, which, as we emphasized before, is not a causal propagator. It is physically 
unacceptable to define a relativistic wave function that would virtually instantaneously spread 
outside the past and future light cone of its region of non-zero support. The Feynman propagator 
is however well suited to do perturbation theory as it is time ordered and hence a Green's function 
(the solution to the Schrodinger equation is a time ordered evolution operator). Despite its obvious 
shortcomings, this is the reason that the Feynman propagator is so widely used when calculating 
Feynman diagrams^. 

We construct real quantum states that include both positive and negative frequency contri- 
butions ibwfc. Our wave function thus becomes the real field itself: 

^{x) = <l,{x) = Mx) + cl,.{x) = I (2^'^)32^^ ['^(^fe' k)e-^->.t+^'^■- + d^i-iOk, k)e+'->^-*+^'^--^ . (3) 



^ There is an additional problem concerning the localization of relativistic particles, i.e.: the standard Newton 
Wigner position operator suflters from non-locality problems that would prevent us from localizing a particle on 
scales smaller than its Compton wavelength. We will not address this problem here but present a solution in [2ll |. 
which further complements a consistent first quantized interpretation of quantum field theory. Additionally, pair 
creation is often argued to fundamentally limit the localization of relativistic particles. This argument does not 
apply however since we are discussing real fields. 
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The positive definite inner product between real fields that includes positive and negative fre- 
quency solutions to the Klein Gordon equation 



22 



probability density: 
1 



(I)\{x)idt4>+{x) - (j)*_{x)idt(j)-{x) 



readily provides us with a positive definite 



(4) 



For real fields, one can easily check that the two currents ([2|) and (j3]) are equal. 

What about the propagation of our quantum states, or of the probability density thus defined? 
Let us recall again that the Feynman propagator propagates complex valued wave functions in a 
manner that is not causal, but it is a Green's function such that it can be used in perturbative 
quantum field theory. There exists, however, a propagator that meets all the necessary require- 
ments: it is real, causal and it is a time ordered Green's function such that we can use it in a 
perturbative expansion as well. This propagator is the time ordered commutator propagator, or 
Wheeler propagator. The Wheeler propagator is the real part of the Feynman propagator, or the 
average of the advanced and the retarded propagator. Only the Wheeler propagator propagates 
the entire real quantum field 4'{x) in a time ordered manner unlike the Feynman propagator. Ul- 
timately, the idea is to formally set up a perturbation theory in quantum field theory using the 
Wheeler propagator instead of the Feynman propagator. 

The structure of the paper is as follows. In section |TI1 we establish our notation and carefully 
develop a consistent first quantization in quantum field theory using real states. Also, we study 
various propagators to examine how they act on quantum states. In section IIII) we consider 
how probabilities in our new setup of quantum field theory are defined. Finally, in section ITVl we 
discuss the role of the Wheeler propagator and see how it causally propagates probability densities. 
Throughout the paper we only consider free, real scalar field theories to keep the conceptual 
subtleties clear and postpone the discussion of interacting field theories to a future publication. 



II. FIRST QUANTIZATION INCLUDING NEGATIVE FREQUENCY 

CONTRIBUTIONS 

First quantization attempts to construct a relativistic quantum field theory for a single particle 
wave function with a probabilistic interpretation just as in non-relativistic quantum mechanics. 
Here, we do not consider second quantization, where creation and annihilation operators are intro- 
duced to conveniently describe multi-particle quantum states. The Lagrangian density for a real 
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scalar field reads: 

£(</), d^4>) = -\ {v'"'0^,4>du4> + ^'0') , (5) 

where we use the Minkowski metric tj^^'^ = diag(— 1, 1,1,1) and where /x is the inverse of the reduced 
Compton wavelength: 

M = = — . (6) 
The Klein Gordon equation of motion follows from equation 1^ as usual as: 

io--^)H^) = o. (7) 

Here, □ is the d'Alembertian. Let us Fourier transform our field as usual as: 



dU 



H^) = J ' (8) 

where we employed the standard notation kx = rj^yk^^x" . Let us set h = c = 1. When we substitute 
equation in the Klein Gordon equation d?]), we find the usual on shell condition (k^)'^ = k'^ +m?, 
such that we can write: 



m= / ^5(fc2 + m2),^(fc)e*'=^ (9) 



We thus see that 0(A;) = 27r(J(fc^ + 'rin?)(f){k). We can clearly see that two separate contributions 
arise if we were to evaluate the delta function explicitly as we have on shell solutions with positive 
and negative frequency: 

/d^A- f d^k 

j^S{k' + m')ei±k^)me"'' = J (2^)32^, '^(^^)"^'^' ' (1°) 

where ook = + m^)^/^ as usual and we defined the contravariant vector: 

k± = {±iOk,k). (11) 

Throughout the paper we will refer to the it indices as positive and negative frequency indices. If 
one insists on keeping a manifestly Lorentz invariant expression, one can leave the delta function 
in equation (jlOp unevaluated. We thus have: 

^{x) = 4>+{x) + <P-{x) ■ (12) 

Since we are studying real fields, we have (j)*{—k) = (j){k) as usual, from which we conclude that 
(p*+ix) = (p-ix). 
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A. Connecting the Schrodinger and Klein Gordon equations 

Although the connection between Dirac's equation and Schrodinger's equation is obvious, one 
should appreciate that similar arguments apply to free scalar fields as well. In non-interacting 
scalar theories, one can factorize the Klein Gordon equation ([7]) as follows: 



i\ + ^A^2_v2j _ ^^2_ V2j = . (13) 

The meaning of the square root of the derivative operator is defined in momentum space. In this 
equation and in the following few we reinstated h and c explicitly again for clarity. Using equation 
([T2D we thus see: 



mc 



indo(j)±{x) = ±mc\ I- — V2 (14) 



These equations are sometimes called "spinless Salpeter" equations. In the non-relativistic limit 
where the relative spatial variations in the field are much smaller than the inverse Compton wave- 
length « (^)^ one obtains: 



±ihdo <P±{x) = m^4>± - ^^^<t>± - ^;^VV± - ■ ■ ■ ■ (15) 



Clearly, the positive and negative frequency parts of the Klein Gordon field satisfy a Schrodinger 
equation and its complex conjugate equation, respectively. We can thus make the following inter- 
pretation: The positive frequency and negative contributions to a quantum field are wave functions 
in the non-relativistic limit. Please note that these considerations only apply to free theories. 

In the literature one can sometimes read confusing statements that the negative frequency 
contributions have negative energy^. The reason for these statements is that the "Hamiltonian" 
appearing in Schrodinger's equation is indeed negative for negative frequency modes. The appro- 
priate Hamiltonian to consider, however, is the Klein Gordon Hamiltonian, which is as usual given 
by: 

H = j d^xdoT^^ = ^Jd^x (^ido(^f + V0 • V0 + m'^cP^^ (16) 
1 f d^k 



2] (2,/('^'' + '"''l'=° 



^ These "negative energy" contributions were considered problematic in the early days of quantum field theory and 
led Dirac to develop his picture of an infinite sea of occupied fermionic states which has now become obsolete. 
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We thus conclude that negative frequency modes have positive energy as they should^. This point 
is of course well known in the literature, but we think it is important to mention it explicitly. 

As (l>+{x) satisfies a Schrodinger equation in the non-relativistic regime and because it has 
positive energies, it was argued some decades ago that (p+ix) perhaps was a promising candidate 
for the relativistic wave function^. It was realized however, that these "wave functions" would not 
spread in a causal manner [l^ . [2^. We take the following point of view: 

• The entire field (pix) and not just its positive frequency part (/>+(x) is a wave function in 
relativistic quantum field theory. The dynamics of the wave function 4'{x) is governed by 
the Klein Gordon equation; 

• The probability density follows from the definition of the inner product on the positive and 
negative Hilbert spaces. 

In other words: we include the negative frequency contributions to a quantum field in what we 
define to be the wave functions in quantum field theory. Note that h explicitly appears in the Klein 
Gordon equation ([7]), as is well known, so we can indeed interpret it as a quantum equation that 
governs the time evolution of the real wave function In the next sections, we properly define 

quantum states, study how they are propagated by various propagators and finally, we associate 
probabilities to the inner product of a quantum state with itself and study how this probability 
density is propagated in time. 



B. Quantum States in Quantum Field Theory 

Let us firstly recall how wave packets are constructed in the conventional approach to quantum 
field theory. Here, the crucial point is that a quantum state is derived from the positive frequency 



contribution to a quantum field only. T 
equation (4.65) of Peskin and Schroeder 



lose states are therefore complex valued. For example, 



331 ] reads: 



(27r)3V2^ 



m\k), (17) 



* We would like to emphasize that we do not introduce fields with negative energy, contrary to the approaches 



taken in e.g. [25l429| . We merely include the Hilbert space of negative frequency solutions to the Klein Gordon 
equation into our description of quantum field theory such that the real field itself can be interpreted as a quantum 
mechanical wave function. 
^ This is one way to view the Feynman-Stiickelberg interpretation of quantum field theory [30l - [33 |. 
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where the subscript "PS" refers to Peskin and Schroeder. In order to better illustrate the point 
we would like to make, we modify this equation slightly to: 



-(j){uJk,k)\ujk,k) . (18) 



Equation ()17p and (jlSp differ in two respects: normalization^ and a notational difference as we 
explicitly include the A;'' dependence. Unlike Peskin and Schroeder, we insist on a Lorentz invariant 
definition for the states. One can express any state in a manifestly Lorentz invariant way as: 

d{LIPS) 4>{k^,k)\k^,k) , (19) 

where the Lorentz invariant phase space (LIPS) element for one particle is in the standard approach 
to quantum field theory given by"^: 

d{LIPS) = -^S{k' + m')0{k') = dk'j^^^Jik^ - u,) . (20) 

It is clear from equations (jlSp and (j20p above that in the Peskin and Schroeder state (|17p only 
positive frequencies +ujk contribute®: 



(21) 



where we added a + sign to the right hand side of this equation to make this explicit. 

Now let us discuss how we define quantum states in our new setup of quantum field theory. 
Given the fact that we can interpret a real scalar field as a wave function in relativistic quantum 
mechanics, we can define a Lorentz invariant state for our field as follows: 

= I -^^(k' + m'mm (22a) 
-^6ik' + m'Wik){k\. (22b) 

Clearly, a quantum state is a superposition of various wave modes |A;) with a certain amplitude 
(f){k) as usual. If one evaluates the delta function the positive and negative frequency modes appear 
again such that the state of a field is a superposition of positive and negative frequency states: 

!</.) = + |<A-), (23) 



^ The Lorentz invariant normalization of quantum states appears for example in [SJl formula 3.4.2. 
^ See for instance Ryder's formula (4.4) [49|. 



The exclusion of negative frequency states is essentially due to Feynman and Stiickelberg [3Ch32| | who understood 
negative frequency states as Hermitian conjugate positive frequency states in the context of charged fields. 
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where: 



r6{k^ + m^)e{±k')m\k) = / Hk±)\k±) , (24) 



and where: 

|fc±) = I iwfc,^). (25) 

It is of fundamental importance to reahze that, unhke the conventional approach taken in the quan- 
tum field theoretical literature as in equation (jl7p . our quantum state (j22p includes the negative 
frequency contributions. This makes our quantum state real. 

The inner product between momentum basis bras and state kets is naturally obtained by de- 
manding that: 

{k±\4>) = (l)±{k) , (26) 
such that: 

{ka\k[) = i27Tf2uJk6^^^{k - k')6ab , (27) 

where a and b are the frequency indices it. This relation was first introduced by Hennaux and 



Teitelboim, who presented a worldline discussion of the super symmetric particle 



221 ]. We can now 



introduce spacetime kets as the covariant on shell Fourier transform of the momentum kets: 

'"^^ = / (0^^^' + m2)e-^-|fc) = + |x_) , (28) 



where: 



\x±) = I ^^5(A;2 +m2)0(±ik°)e-*'=nfc) = / . .„ e-'''^'-'\k±) (29a) 



(27r)3 ' ^ \ ^ ' ' J (27r)32wfc 

vm^' '•'^^ ■ 

The sign in both exponents is chosen to give the standard sign for the Fourier transform in equation 
([8]) as we show below in equation (|34p . Equation ()28p merits another remark. If we compare this 
equation for \x) with equation (|22ap for we can identify 4>{k) <^ e~^^^. In other words: the state 
\x) corresponds to a superposition of plane wave modes, i.e.: a non-normalizable "wavepacket" of 
plane waves. The inverse Fourier relations are given by: 

|A;±) = T J d^xe^^i^i^olxi) (30a) 

(A;±| = ± j d^xe-'''^^ido{x±\ . (30b) 
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Here, we defined a time derivative that acts in the following way on two test functions f{x) and 
g{x): 

/(x) dog{x) = f{x)dog{x) - g{x)dof{x) . (31) 
Having carefully defined the spacetime states, we find: 

cj){x) = Mx) + = {x+\<i)) + (x_|</.) = (x|0) . (32) 

We can furthermore derive: 

\(^±) = ^ j d?x(^±{x)ido\x±). (33) 

Using (j29p and the positive definite momentum space inner product (j27p one easily derives: 

{x±\k±) = e^'^i^ (34a) 

{k±\x±) = e-^*^±^ (34b) 

{x±\k^)=Q (34c) 

{k±\x^) = Q. (34d) 

C. Inner Products and the Identity Operator 

We can also compute various position space inner products. The inner product between the 
positive and negative frequency modes separately yields the two Wightman functions: 

A±(x-x') = ±i(x±|x'±) = ±i I -^,5(fc2+m2)0(±A:°)e^'^("-"') = ±i f ^ik±(x-x') _ (35) 

J (2vr)'^ J (27r)'^2a;fc 

The ibi in the definition of the Wightman functions becomes apparent later when we discuss 
propagators. Moreover, it is important to stress that our it frequency indices are very different 
from the Schwinger-Keldysh indices often used in the in-in formalism in out-of-equilibrium quantum 
field theory. We also find: 

{x±\x'^) = Q. (36) 

The Wightman functions furthermore obey: 

A+ (x - x') = - A" (x - x) (37a) 
[A±(x-x')]* = A^(x-x'). (37b) 
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The positive definite inner product between two arbitrary states follows from equation (j22p as: 

^6{k^ + m^)x*{k)m = j (2^)32^J (x|fc+)(fc+|0) + (x|fc-)(fc-|0)] (38) 

From the equation above we read off a Lorentz covariant identity operator: 

i = f ^5{e + m')\k){k\=U + L, (39) 

where: 

which we can easily verify by for example applying it to \k±). To derive the position space identity 
operator let us consider objects like {x±\x'j^ derived in equation (j35p . Inserting the following 
operator leaves these functions unchanged: 

4 = ±y d?x\x±)idQ{x±\ , (41) 

such that the position space identity operator on the total space of solutions is given by: 

d^x\x+)ido{x+\ - \x-)ido{x_\ . (42) 

D. Propagators, Propagation Laws and Composition Laws 

In this subsection, we derive various propagators for our quantum field. Intuitively, it is im- 
portant to realize that propagators only propagate a state but do not contain information about 
the state itself. In other words: the wave mode \k) of a certain state has an amplitude (j){k) as 
in equation ()22ap . Wave packets can thus be constructed by the principle of superposition. Prop- 
agators do not contain any of this information. The amplitude of a certain wave mode \k) for a 
propagator is just the plane wave (t){k) = e"**^^. A propagation law captures how a propagator acts 
on a certain wave function, whereas a composition law captures how a propagator acts on another 
propagator. 

1. Wightman propagators 

We have already encountered the two Wightman functions in equation ()35p . given by: 
A^(x — x') = ibi(x±|x'_j_). What we discuss here is how they act on the wave functions. By making 
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use of our identity operator (j4ip one can straightforwardly show that: 

(p±ix) = ±i j d^x'{x±\x'^)dt'(l)±ix') = j d3x'A±(2; - x')dt'<l)±{x') . (43) 

The propagators for positive and negative frequency wave functions are the positive and negative 
frequency Wightman propagators, respectively. Clearly, a wave function at a time t' is transformed 
to a wave function at another time t, where at the moment t can be either earlier or later than 
t' . It is important to realize that information about the state, i.e.: the specific superposition of 
wave modes in (|22p. is always contained in (j)±^ whereas the propagation of a state is governed 
by the Wightman functions. Wightman propagators can also act on the complex conjugate wave 
functions: 

(/.^(x) = d^x'(i)*^{x')dt'{x'^\x±) = j dV</.i(x')5t/A^(x' - x) . (44) 
Here, we have made use of equation (|35|) . Moreover, we can easily see that: 

d^x'A±(x - x')dt'<j)^{x') = (45) 
The Wightman propagators satisfy the following composition laws: 

dVA±(2; - x")dt',/\'^{x" - x') = A±(3; - x') , (46) 



which we prove straightforwardly by making use of the integral representation (j35p . It will come 
as no surprise that the mixed composition laws for the Wightman propagators vanish: 

d^x" A^ (x - x") dt" A^ (x" - x') = . (47) 

2. Feyman propagator 

The Feynman propagator ensures that the "right to left" direction of propagation in equation 
3|) coincides with the direction of increasing time: 

9{t - t')(l)+{x) = j d^x'Apix - x')df(j)+{x') (48a) 

e{t' - t)(l)l{x) = j d^x'(l)l{x')dt'AF{x - x') , (48b) 

where the Feynman propagator is the time ordered positive frequency Wightman propagator: 

A^(x - x') = iT[{x+\x'^)] = ie{t - t'){x+\x'+) + ie{t' - t){x'+\x+) (49) 
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Here, T denotes the time ordering symbol. This form shows exphcitly that the Feynman propagator 
describes the on shell time ordered propagation of the positive frequency modes of a field. In other 
words, the Feynman propagator is a retarded propagator for positive frequency modes. Also, we 
see that the Feynman propagator is an advanced propagator for the complex conjugates of positive 
frequency modes. Since we know that (^4.(2;) = (j)-{x), we conclude from equation (|48bp that the 
Feynman propagator is an advanced propagator of the negative frequency contributions too. We 
can bring equation (|49|) in the conventional form by also Fourier transforming the Heaviside step 
functions: 

dE f d^k 



I ^-iE{t-t')+ik+(x-x') _ ^ ^-iE(t-t')-ik+{x-x') 



27r J {27rf2L0k [E + ie" E - ie 

d^k 1 



^ik{x-x') ^ (-^Q) 



(27r)4 A;2 + m2 - ie 

where we shifted the E integration appropriately in the second line. Alternatively, we can rewrite 
the equation above as: 



Air(x-x') = P 



d'^k 1 



^ik{x—x') 



1 ^A^'+^'y'^^-^'^^ (51) 



(27r)4 fc2 + m? 

where we have made use of the Sokhotski-Plemelj or Sokhotsky-Weierstrass formula, also known 
as Dirac's rule: 

^ P- T ivr5(^) , (52) 



where P denotes the Cauchy principal value. The anti-Feynman or anti-time ordered propagator 
(T) is now given by: 

Aaf{x - x') = iT[{x+\x'_^)] = ie{t' - t){x+\x'^) + ie{t - t'){x'^\x+) (53) 

^ f 1 ik(x-x') 

J (2^)4 it2 + m2 + ie 

Let us now discuss the various composition laws for the Feynman propagator. Firstly, a Feynman 
propagator can act on another Feynman propagator: 

J dVAi.(x-x")S"Ai.(x"-x') = 0(i-i")^(i"-i')A+(x-x')+e(t"-t)0(t'-t")A-(x-x'), (54) 

where it is crucially important to make use of the integral representation of the Feynman propagator 
in equation ()50p in order to take the distributional character of the propagator properly into 
account. If we time order the expression above, i.e.: we require the latest to the left convention, 
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we find the composition law for the Feynman propagator: 



T 



d^x"AF{x - x")dt"AF{x" - x') 



(55) 



e{t - t') / <rx"AF{x - x")dt.AF{x" - x') + e{t' - t) / d^x"AF{x' - x")dt"AF{x" - X 



Af{x - x') j„ + 1:Af{x - x') ^„,^ + \^f{x - x') 



t"=t' 



t"=t 



Also, we can study how a Feynman propagator or anti-Feynman propagator acts on our Wightman 
functions. By making use of equation ()50p and of several complex contour integrations one can 
show: 



(56a) 
(56b) 



d^x"AFix - x")dt"A^{x" - x') = ±e{±t ^ t")A^{x - x') 
J d^x"AAF{x - x")S"A±(x" - x') = TO{Tt ± t")A^{x - x') . 



Apart from being a propagator, the Feynman propagator is also a Green's function for the Klein 
Gordon equation. By making use of equation (|5Up one can verify that: 



(-a^ + m'^)AF{x - x') = <5(^)(x - x') , 
such that the solution of (—9^ + m^)(/>+(x) = J{x) is: 



(57) 



(p+{x) = /(x) + / d^x'AFix - x')J{x') . 



(58) 

Here, (f)^{x) is a complex valued solution to the homogeneous Klein Gordon equation, and J is 
a complex valued source. Note that the Feyman propagator is complex valued and is therefore 
certainly not the Green's function of the real scalar field (j). It can however be the Green's function 
of e.g. the positive frequency contribution to our real field. 



3. Commutator Propagator 

So far, we have introduced the two Wightman and Feynman propagators. The former propagate 
the positive and negative frequency contributions to a real field, whereas the latter act as retarded 
and advanced propagators for these contributions, respectively. There is, however, another prop- 
agator worth discussing. It is the commutator propagator and it propagates the entire field, i.e.: 
both its positive and negative frequency contributions: 

(p{x)= f d^x'Ac{x-x')dt'(p{x') (59) 
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where the commutator propagator Ac{x — x') is given by: 

Ac(x — x') = irf^ {xa\x'ij) = i{xj^\x'j^) — i{x^\x'_) = A^(x — x) + A~(x — x') . 



Its integral representation reads: 



Ac(x — x') = i 



(60) 



(61) 



where e(A;'') is sign function as before. It is important to reahze that the commutator propagator 
satisfies the following relations: 



/\c{x-x\^^, =0 
dMx - x\^^, = 5\x - x') . 

These identities are important for the following reason: 

- / d'^x'(l){x')dt'Ac{x-x') 

t=v J 

{x,t'). 



(62a) 
(62b) 



lim(j){x)= / d x'Ac{x — x')dti(l){x') 



t=t' 



(63) 



Thus at time coincidence no propagation of the field has occurred. The commutator propagator 
satisfies the following composition law: 

Acix -x')= [ d^x'Acix - x") dt"Acix" - x') . (64) 



4- Wheeler Propagator 

Clearly, the commutator propagator is not a Green's function as it does not satisfy the inho- 
mogeneous Klein Gordon equation. Let us therefore define the Wheeler propagator as the time 
ordered commutator propagator: 



Aw{x - x') = \t[Ac{x - x')] = 'J^^Acix - x') 



(65) 



The last identity follows from the anti-symmetr y o f the commutator. The Wheeler propagator 
is thus symmetric in its arguments^. We follow [37|] in calling the time ordered commutator the 



Wheeler propagator^ as it appears implicitly in Wheeler and Feynman's absorber theory 



Wheeler's two point function is also the average of the advanced and retarded propagator 
1 



39 



38|. 
4qI: 



Awix — x') 



{Aa{x - x') + Ar[x - x')) 



(66) 



9 



This propagator appears in a classical context in [35], |3£ 



We do not agree however with the statement made in [33 that the Wheeler propagator lacks on shell propagation, 
as is clear from equation (|69p . 
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where: 



Ar{x - x') = 9{t - t')Acix - x') 
Aa{x - x') = -e{t' - t)Ac{x - x') . 

Moreover, the Wheeler propagator satisfies the fohowing propagation law: 
e{t - t') 



d-^x'Awix-x')dt'(l){x'). 



(67a) 
(67b) 

(68) 



Hence, like the commutator propagator, the Wheeler propagator propagates the entire real quan- 
tum field. Note that this is natural since Wheeler's two point function is real unlike the Feynman 
propagator. The Wheeler propagator has the following integral representations: 



Aw fx — x') 



2.1 (27r)3 



J,k{x—x') 



2 J (27r)3 

1 (■ d^fc 

2 J (2^ 



1 1 

k"^ + m? + ie k"^ + m? — ie 



J,k(x—x') 



(69) 



where on the second line we used the Fourier representation of the Heaviside step function. The 
Wheeler propagator is thus given by the real part of the Feynman propagator: 



Aw{x -x') = ]^ {Af{x - x') + Aaf{x - x')) = P 
= ^Af{x-x'). 



d'^k 



1 



(27r)4 A;2 + 



^ik{x—x') 



Here, we used Dirac's rule in equation ()52p again. Clearly, this would equally well al 



(70) 



ow us to refer 



4l[ |. Let us now 



to the Wheeler propagator as the "principal part propagator", see for instance 
discuss the composition law for the Wheeler propagator. A Wheeler propagator acts on another 
Wheeler propagator as follows: 



1 



d^x"Aw{x - x")dt"Aw{x" - x') = -Ac{x - x')e{t - t")e{t" - t') . 



(71) 



Here, we made frequent use of equation ()56p . Like the composition law of the Feynman propagator, 
we should time order the expression above. We thus derive: 

1 



T 



d-^x"Aw{x - x")dt"Aw{x" - x'] 



-Aw{x-x')e{t-t")e{t" -t'). 



We can recognize this equation as a composition law if we require t < t" < t' U t' < t" < t: 



Aw(x-x') = 2T 



d^x"Aw{x - x")dt"Aw{x" - x') 



(72) 



(73) 



t<t"<t' u t'<t"<t 
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The Wheeler propagator is not only a propagator, it is also a Green's function of the Klein Gordon 
equation because of the time ordering. One can again easily see that: 

{-d^ + m'^)Aw{x - x') = - x') , (74) 

where we have made use of equation ([69]) . The solution of (—9^ + m'^)(j){x) = J{x) is thus: 

<p{x) = 4P{x) + j dVAvy(x - x')J{x') . (75) 

Here, 4P{x) is a real solution to the homogeneous Klein Gordon equation. Finally, note that the 
Wheeler propagator is the Green's function of the real scalar field cj)^ and not just of the complex 
valued positive frequency contribution as the Feyman propagator is. 



III. PROBABILITY DENSITIES IN QUANTUM FIELD THEORY 

Given our wave function interpretation of the real quantum field (p as presented above, it is 
natural to define probability densities analogous to quantum mechanics. In quantum mechanics, the 
absolute value squared of a quantum mechanical wave function represents a probability distribution 
function of a single quantum mechanical particle, whose dynamics is governed by Schrodinger's 
equation. The absolute value squared of a quantum mechanical wave function can be used as a 
probability measure because firstly it is positive and, secondly, it can be normalized to unity when 
integrated over its domain. The natural candidate for a probability density in relativistic quantum 
mechanics stems from the inner product between a quantum state with itself as we show shortly. 
The probability densities thus obtained are positive and can be normalized to unity too. 

In the conventional approach to quantum field theory this natural generalization of quantum 
mechanics is problematic as only the positive frequency modes are included in wave packets and 
the negative frequency modes do not contribute. The propagation of these wave packets is, as 
we have seen above, governed by the Feynman propagator. The wave packets therefore spread 
non-causally. We solve this problem by properly including the negative frequency contributions in 
the wave packets such that their propagation is mediated by Wheeler's propagator. We therefore 
have a normalized probability density at our disposal that spreads in a causal manner. This allows 
us to build a consistent first quantized picture in quantum field theory that naturally generalizes 
quantum mechanics to a relativistic setting. 
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A. A Normalized Probability Density and the Hadamard Propagator 



In our setup of quantum field theory we require: 
Inserting the identity operator from equation ()42p yields: 



1 /• „ ^ ^ 

- / d-'x{(t)\x+)idQ{x+\4>) - {(j3\x^)idQ{x^\i 

- / d-^ x(l)\{x)ido())+{x) - (l)*_{x)ido(l)-{x) 



(76) 



(77) 



d^k 



2 J (27r)32wfc 

1 r d^fc 



6ik'^ + m^)\(j){k)\'^ . 



(27r)3 

The momentum space realization of this inner product is discussed for example in 



42l | where 



231 



the author considers a stochastic interpretation of the fields. Halliwell and Ortiz discuss in 
both the position space and momentum space realization in the context of the sum over histories 
interpretation of quantum field theory. Their work was inspired by earlier work of Henneaux and 
Teitelboim 



22( 1 . In a slightly different form the inner product can also be found in [2j], also see 



43l |. The factor of 1/2 that appears in ()77p is motivated by the canonical form of the Hamiltonian 
in equation (jl6p . It is clear from equation (j77p that plane waves cannot be normalized. One really 
needs a wave packet, i.e.: a superposition of various wave modes with a certain amplitude, such 
that the integrals in ()77p converge. Note that a single plane wave does not meet this requirement 
and is therefore not a suitable wave packet. We can thus identify a probability current: 



(t)\{x)id ^(t)+{x) - (t)*_{x)id ^(t)^{x) 



which is conserved for free theories in the usual sense: 



5^j;(x) = o. 



(78) 



(79) 



where we have used the Klein Gordon equation. We identify the 0th component of the probability 
current as the probability density, given a wave function (j){x): 



pp{x) = J^{x) 



1 



(i)*Ax)idt(l)+{x) - (I)*^{x)idt4>-{x) 



(80) 



From the third line in equation ()77p one can clearly see that the probability density is both real 
and positive. Given a spacelike region A, the probability of finding a particle in this region for 
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some observer thus follows as: 



1 



<l)l{x)idt(l)+ix) - (j)*_{x)idt(l)-{x) 



(81) 



Y'xJ^ix) = j d^xpp{x) = 2 y 
A A A 

We already tentatively used the word particle in defining the probability above. We make this 
statement more precise shortly. Let us however first define the statistical propagator, or Hadamard 
propagator as: 



A/f(x — x') = zA^(x — x') — iA {x — x') = — (x+|x'_,„) — 



(82) 



5(fc2 + m2)e^*^(^-^') 



^ik+{x—x') _|_ ^ik—{x—x') 



(27r)3 ' ' J (27r)32wfc 

We can straightforwardly relate the Hamadard propagator to the imaginary part of the Feynman 
propagator: 



Ah{x - x') = -2QAf{x - x') . 



(83) 



The Hadamard propagator satisfies the following composition law in terms of the commutator 
propagator: 



/ 



d^x'Acix - x")dt»AH{x" - x') 



Affix - x') . 



(84) 



The reason for introducing the Hadamard propagator now becomes apparent. We can now show 
that: 



^{q^lcP) = ^ I d^fdV cPix)dtAH{x - x')dt'cl){x') 



1 



(85) 



t=t' 



This inner product can also be found in 



23l |. The Hadamard propagator thus offers an alternative 
way of evaluating the probability for finding a particle in our spacelike region A for a particular 
observer: 



1 f ^-^ ^-^ 

P=- d^xd^x (j){x)dtAH{x - x')dt'(l){x') 



(86) 



t=t' 



The probability density J(f (x) is not Lorentz invariant, which is especially clear from the position 
space expression in equation (jSOp . The probability density current Jj^{x) is, however, Lorentz 
covariant. Also, the normalization of our probability to one, where J^(x) is integrated over its 
entire domain, is Lorentz invariant. The reason is that the inner product integral contains a Lorentz 
covariant vector valued volume element, which contracts with the current to give a Lorentz invariant 
result. 
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The probability current J^{x) in equation (j78p merits a final remark. For real fields, we have 
that 4>~{x) = (j)*^{x). Clearly, the probability current is in that case equal to the ordinary Klein 
Gordon current, given by: 



iv^^(x) = </.;(x)z9^0+(x) 



(87) 



where the negative frequency contributions in (|78|) are mapped to the positive frequency Hilbert 
space. Therefore, equation ()87p is a probability density too. Such an interpretation can for example 



be found in 



13, 



18|, 



34( 1 . Although J^(x) and N^^ (x) are equal, there are profound differences 
between the dynamics of the wave function (j){x), which is causal, and of (j)-^.{x), which is not causal. 
Therefore, the positive frequency contribution to the real field on its own cannot be a relativistic 
wave function. 



B. Particle Interpretation 

In non-relativistic quantum mechanics, a wave function describes the probability density of a 
single particle. Motivated by the discussion above, it is natural to interpret the probability density 
in equation (jSOp as the number density for a single particle too: 



N = \m) = \j d3f0;(x)iao</.+ (x) -r_(^)i5o</.-(x) = 1. (88) 

In other words, it is the superposition of positive and negative frequency wave packets, normalized 
to unity, that defines one particle. Hence, our particle interpretation is in agreement with the 
particle interpretation as advocated in non-relativistic quantum mechanics. Take for example the 
Gaussian ground state of a simple harmonic oscillator. Clearly, it consists of the superposition of 
many wave modes but is nevertheless generally regarded as a single particle quantum mechanical 
wave function. This differs from the particle interpretation that is usually adhered to in quantum 
field theory, where one mode \k) = aj^|0) is interpreted as one particle, see e.g. j^l- It also differs 
from various claims made in the literature that the wave function for a single relativistic particle 



does not exist in quantum field theory 



11 



19|, 



201, 



44 



46|. 



The particle number in the sense defined above in equation ()88p is only conserved in free theories. 
As for the electric charge density, and for the probability density, we can define a Lorentz covariant 



We are able to adhere to a single particle wave function interpretation since in our framework wave functions 
propagate causally. As mentioned before however, there are still problems related to non-locality in the standard 
treatments of position operators. This issue will be resolved in [211 ]. 
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particle density current: 

Nf. = ^ct>li^)id^M^) - l(t>*.ix)id^^.{x) . (89) 

For complex fields for instance the charge current and the energy momentum current are conserved 
because the Lagrangian following from ([5]) is invariant under U{1) and it is translationary invariant. 
The particle number density current defined in equation (|89|) is not protected by such a symmetry 
in interacting theories. 



IV. CAUSAL PROPAGATION OF QUANTUM FIELDS 
A. The Hadamard Composition Law and the Wheeler Propagator 

In the previous sections, we have developed a first quantized picture for free scalar quantum 
field theories. In this final section, we rigorously show that probabilities, defined in e.g. equation 
(j86p where we made use of the Hadamard propagator, propagate causally. 

Let us firstly remind the reader again that the Feynman propagator is not a causal propagator. 
Suppose that, at some initial time t', one specifies a certain state defined in equation (|22ap . 
Let us now recall equation (j48ap : 

e{t- t')(j)+{x) = J d^x'Apix - x')dt'ct)+{x') . 

The Feynman propagator has a non-vanishing support outside the past and future light cone 
of x' and hence we see that it propagates the positive frequency contribution to our real field 
(f> instantaneously outside the future light cone of its region of non-zero support. The Wheeler 
propagator stems from the commutator, and thus only has non- vanishing support inside the past 
and future light cone of x' . If we recall equation ([68]) : 

e{t - t') 



2 



x) = / <rx'Aw{x - x')dt'(t>{x') , 



we conclude that (t){x) will only be non-zero if x is the past or future light cone of the region of 
non-zero support of 4>{x'). Qualitatively, this is depicted in figure [TJ 

The probability for some observer for finding a particle in some spacelike region can be 
calculated by making use of the Hadamard propagator, as we recall from equation ([86]) : 



1 r ^-^ -f-^ 

P=- d^xd^x' <t){x)dtAH{x - x')dt'^{x' 

A 



t=t' 
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Figure 1: Qualitative picture of the causal propagation of a quantum field mediated by the Wheeler 
propagator. Note that the evolution ensures that no contributions propagate beyond the causal region of 
the initial state. In other words, outside the causal future of the wave packet there is perfect destructive 
interference between the positive and negative frequency modes. 

We can derive the following composition law for the Hadamard propagator: 

--^e{t-ty)e{t' -ty,)i^H{x-x') = j d?y'd?f^w{x-y)dt^l^H{y-v')\Awy-x'). (90) 

Here, we made use of equations (|56|) and (|7U|) . We have thus expressed the Hadamard composition 
law in terms of the Wheeler propagator rather than the commutator propagator as in equation 
()84p . Also, we defined y = {yt,y). Alternatively, we can absorb the minus sign in front of equation 
(PU|) in the composition law as follows: 

^e{t-ty)e{t' -ty,)AH{x-x') = J d^y'd^f Aw{x - y)dt^AH{y - y')dt^,AAw{y' - x') , (91) 

where the anti- Wheeler propagator is defined as the anti-time ordered commutator propagator: 

Aaw{x - x') = ^TAc{x - x') = -Aw{x - x') . (92) 
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If both t and t' are later than ty and ty' , the composition law (j9ip takes a rather neat form: 



{t,t'}>ty ■ (9^) 



{t,t'}>ty, 

We already showed before that the Wheeler propagator governs the causal dynamics of the entire 
wave function (j){x). Now, together with (|86p . equation (|93p shows explicitly that probabilities 
propagate causally too as the time evolution of the Hadamard propagator can be expressed solely 
in terms of the causal Wheeler propagator. A final technical remark is in order. The factor of 
4 that appears on the right hand side in equation (|93p arises from our definition of the Wheeler 
propagator as a factor of half times the time ordered commutator ()65p . such that the Wheeler 
propagator is a Green's function. If we would not have inserted the factor of 1 /2 in our definition, 
the Hadamard propagator would satisfy a neater looking composition law. 

B. Example I: a Localized Wave Function 

In order to illustrate the discussion above, let us study the propagation of two wave packets: a 
specific (p+ix) and <^(x). We consider the following wave function (p: 

^(^) = f + gifc-(x-x')l _ (94) 



(2^)3 2 

We have thus chosen the following wave packet: 

4>{k±) = Uke-'"^^' , (95) 
The second wave packet we wish to consider is (j)+{x) which can be read off from equation ()94p : 

These wave packets do not give rise to a probability density that can be normalized according to 
equation ()76p . The reason for introducing these wave packets nevertheless is that they are perfectly 
localized initially: 

Hx)\t=t' = ^^^Hx - x') (97a) 

Mx)\t=t' = ls^'H^-^')- (97b) 

These (extremely well) localized wave packets allow for a clean comparison between the form of 
the wave packets for t > t'. The conventional Newton- Wigner wave packet ^] is not given by 
equation but rather by: 

^f^ik) = a;|e-^'=+^' , (98) 
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and, consequently, gives rise to a wave function that is not perfectly localized initially. In fact, the 
position space expression for the Newton- Wigner wave packet yields at t = t': 

<Ar(^)|*-t' = 7 3^ f ' KM\x-x'\\) , (99) 

247r2r(i) \\\x-x'\\J 4 

where Ku{x) is a modified Bessel function of order v. This expression reveals that the Newton- 

Wigner wave packet is not perfectly localized initially. 

All we need to do next, is perform the Fourier integrals in (|94p and (|96|) . We make use of the 

form of the two Wightman propagators (j35p and introduce two e pole prescriptions to regulate the 

integral. Let us for simplicity consider a massless field such that Uk = \\k\\. Our wave function is 

given by: 

^{x) = Cdf, [A+(x - x') + A-{x - x')] = dt [A+(x - x') + A-{x- x')] (100) 



A 

47r2 
ijt - t') 
27r2 



-{t -t' - ie)2 + ||f - x'P -{t -t' + ieY + ||f - f 



1 1 



. [-{t -t' - k)2 + \\x - x'W^f [-{t -t' + ze)2 + \\x- x'\\^f_ 
In the first equality we have chosen our state such that is normal to the spatial hypersurface. 

Moreover, the equations on the second and third line are to be interpreted in the distributional 

sense where e — t- 0. Most importantly, our wave packet spreads in time and gives a non-vanishing 

contribution only on the light cone. This is to be contrasted with the result for (j)+{x) which results 

from propagation forward in time with the Feynman propagator: 

^+{x) = dtA+{x-x') = lit^^U (101) 

27r2[-(t-t'-ie)2 + ||f-x'||2]2 ^ ^ 

This wave packet has non-vanishing contributions outside the light cone too. We conclude that 
although initially both wave packets are extremely well localized as delta functions, at later times 
spreading of the wave packets has occurred such that our wave function (j){x) has a non-vanishing 
contribution only on the light cone whereas 4>+{x) has not propagated in such a causal manner and 
has a non-vanishing contribution outside the light cone too. This example thus shows that only 
<j){x) can be a wave function consistent with the relativistic principle of causality. 

C. Example II: a Localized Probability Density 

We can also consider relativistic wave packets that can be properly normalized according to 
equation ()76p . Let us now consider another wave function (j){x) which is given by: 

(l)(x) = A /" ^J/" Lik+ix-x'-iO ^^ik-ix-x'+ii)] _ /^02) 

V-Af J (2vr)'^2a;A; L J 
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Note that with this normahzation 4>{x) and are Lorentz scalars. This time, we have thus chosen 
the following wave packet [48[: 



^{k±) = ^e-*^±(-''±*«) , (103) 



where ^ is a timelike vector with > ||^|| such that the momentum space integral converges 
and where M is a normalization constant. For simplicity we will consider = (e, 0, 0, 0) only. 
In particular, we will be interested in the limit where e is small. In principle, however, we can 
generalize this analysis to more general states where is arbitrarily large and where ||^|| 7^ as 
long as the convergence condition > \\^\\ is satisfied. Physically, states with ||,^|| 7^ have a 
non-zero initial velocity of the wave packet. Also, determines the width of the wave packet. It 
is interesting to note that for m ^ 0, x' = and for small ||A;|| the wave packet in eq uation ()103p is 
a Gaussian. For a more general discussion of the dynamics of wave packets, see 2l|]. Let us again 
specialize to massless scalar fields for which ojk = \\k\\- We can normalize our state by making use 
of equation ([76]) : 

We can evaluate the Fourier integrals that appear in equation (jl02p to find our wave function in 
position space: 

^^""^ " TT[-{t-t' -ie)^ + \\x-x'\\'^] ^ TT[-{t-t' + ie)^ + \\x-x'\\^] ■ ^^^^^ 



We are interested in the probability density which we defined earlier in equation ()80|) : 
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, , 4e3 [{t-t'y + \\x-xT + e , 

Pp(x) = ^ K . (106) 



7r2 [-{t -t' - ie)2 + ||f - x'WY [-{t -t' + ie? + P 



As a simple computational check of the normalization of our probability density, one can easily 
verify that: 

^xppix) = 1 , (107) 



for arbitrary times as it should. In figure [2] we show the real wave function in equation (|105p at two 
different times. In ordinary quantum mechanics, one cannot so easily visualize the wave function 
as it is complex valued. In figure [3l we show the resulting probability density in equation (|106p . 
Note that it is a spherical wave travelling at the speed of light. The reason for considering a wave 
function of the form (|102p is that initially its corresponding probability distribution is extremely 
well localized: 
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Figure 2: The wave function given in egua- 
tion il05\) at t/eQ ~ 1 and t/eg ~ 2. In non- 
relativistic quantum meehanies, it is hard to 
visualize the wave function as it is complex val- 
ued. Here, however, our wave function is the 
real Klein Gordon field. We used e/eo = 0.1. 



Figure 3: The probability density in equation 
\10b]) resulting from the wave function depicted 
in figure\^ Note that we multiplied by the sur- 
face of the spherical shell. We used e/eo = 0.1. 



Sending e to zero yields: 



lim pp{x) 



1 



e->0 ' ' ' 2'k\\x — X'P' 

Interestingly, the probability density for t — t' > is given by: 

1 



(109) 



limpp(x) 



AttWx — T 



•/||2 



[6 {\\x - x'W -it- t')) + 6 {\\x- x'W + {t- t'))] . 



(110) 



For massless scalar fields, we thus observe that the probability density propagates exactly on the 
light cone. Furthermore, it is important to realize that although 4>+{x) evolves with the Feynman 
propagator and consequently spreads outside the light cone of its initial region of non-zero support, 
the probability density at late times behaves perfectly causal. The reason is that our probability 
density given in equation ()80p and the Klein Gordon current in equation ()87p are equal. This 
is of course not surprising. The reason for introducing the Wheeler propagator and real valued 
relativistic wave functions has not been to prove causality at the level of probabilities as expectation 
values in quantum field theory are already perfectly causal. The main point of our formalism, 
however, is that in-out amplitudes themselves can be treated causally too. Finally, we would like 
to stress that our single particle real valued relativistic wave function in equation (jl02p is localized 
and can be properly normalized. This is to be contrasted with the states constructed by Newton 
and Wigner 47| as the eigenstates of the Newton- Wigner position operator. 
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V. CONCLUSION 

We have carefully developed a first quantized interpretation for real scalar fields in non- 
interacting scalar field theory: 

• The wave function of quantum field theory is the entire real field (j){x), which generalizes the 
wave function of non-relativistic quantum mechanics; 

• The dynamics of the wave function is governed by the Klein Gordon equation, which gener- 
alizes the Schrodinger equation in quantum mechanics; 

• The probability density follows from the innerproduct on the Hilbert space of positive and 
negative frequency contributions to (j){x). The probability density is of course normalized to 
unity when integrated over its domain of non-zero support. 

By making use of the stress-energy tensor for a real field we also showed that the negative frequency 
solutions have positive energy. The particle interpretation made in a first quantized picture of 
quantum field theory differs from the particle interpretation that is usually advocated in quantum 
field theory. We interpret a wave function normalized to unity as the quantum mechanical wave 
function of a single particle. Due to the superposition principle, this wave packet thus contains 
many wave modes. In conventional quantum field theory, one wave mode is usually interpreted 
as one particle, see for instance [4^ . The dynamics of the single particle wave function (j){x) is 
governed by the Wheeler propagator, which is the time ordered commutator propagator. The 
Wheeler propagator has several advantages when compared to the Feynman propagator: 

• Unlike the Feynman progagator, the Wheeler propagator is causal. It does not have any 
contributions outside of the past and future light cones of its arguments x and x'; 

• The Wheeler propagator is real and propagates the entire real field (j){x), while the Feynman 
propagator is complex valued and only propagates the positive frequency contribution (j)^{x); 

• Both the Wheeler and the Feyman propagator are Green's functions, which makes them par- 
ticularly convenient to develop perturbation theory. However, only the Wheeler propagator 
is the Green's function of the entire real field <j). 

Precisely because the Wheeler propagator or the commutator propagator governs the causal prop- 
agation of the real field (j){x), we can interpret the field itself as the wave function in relativistic 
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quantum mechanics. This is to be contrasted with the function (/>+(x) that is virtually instanta- 
neously propagated outside the past or future light cone of its region of non-zero support by the 
Feynman propagator. 

In our approach the probabilities are calculated using the Hadamard propagator. We only need 
the equal time Hadamard propagator, so it plays no role in the dynamics of our scalar field. The 
Hadamard propagator satisfies a composition law that can be expressed solely in terms of the 
Wheeler propagator. We have thus shown that probabilities in our form of the in-out formalism 
propagate causally. In other words: causality is restored in the in-out formalism of quantum field 
theory by properly including the negative frequency solutions to the Klein Gordon equation. 

We would also like to mention again that expectation values of observables in the first quantized 
picture of quantum field theory will not be affected by whether one makes use of the Feynman 
propagator or of the Wheeler propagator to calculate them. So, standard quantum field theory 
in the in-out formalism is causal when expectation values of observables are considered, but not 
when amplitudes are calculated. However, by making use of the Wheeler propagator, we have 
been able to implement causality explicitly at the level of the dynamics of the fields, unlike in the 
conventional picture of first quantized quantum field theory. 

The first quantized picture of quantum field theory we discussed here is only complete if we can 
derive a local position operator. This issue will be addressed in a separate publication [21[. As said 



at various places throughout the paper, we have only considered free real scalar fields. Of course, 
we are in the process of generalizing our setup to interacting scalar fields and set up a perturbative 
expansion using the Wheeler propagator to calculate expectation values. Also, it would be very 
interesting to generalize our setup to e.g. de Sitter spacetime. 
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